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1*  Introduction:  In  a previous  report  [1]  an  exhaustive  discussion  of  the  qualitative 
behavior  of  the  solutions  to  the  equation 


- B(x) 


led  to  a complete  classification  of  stable  and  unstable  characteristics.  It  was 
shown  that  system  (1)  can  never  be  asymptotically  stable  and  that  a necessary  and 
sufficient  condition  for  the  origin  to  be  stable  is  that  B(x)  be  of  the  form 


B(x)  - c xDx 


where  the  matrix  D has  complex  eigenvalues.  In  this  report  we  consider  a more 
general  differential  equation  of  the  form 


x - Ax  + B(x) 


with  both  linear  and  quadratic  terms.  While  the  qualitative  behavior  of  system  (2) 
was  exhausted  [1]  using  a few  simple  prototypes,  the  addition  of  the  linear  part  in 
(3)  results  in  much  more  varied  behavior:  a simple  classification  will  not  be 
complete;  a complete  classification  will  not  be  simple.  Our  main  Interest  is 
consequently  in  the  asymptotic  stability  in  the  large  of  equation  (3).  The  princi- 
pal result  presented  here  is  that  the  conditions  — [ — ! R ““ 1 

(1)  A is  a stable  matrix  jj  J ^ i fr;  •.[ 

(ii)  B(x)  - cTxDx  2 t fc 


% 
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(111)  the  matrix  D has  complex  eigenvalues 
(lv)  the  matrix  A has  complex  eigenvalues 

(v)  the  matrix  A + Xd  has  the  eigenvector  for  some  real  X and  Its 
corresponding  eigenvalue,  y,  Is  negative 

are  necessary  and  sufficient  for  the  asymptotic  stability  in  the  large  of  the  equili- 
brium state  x ■ 0.  In  section  2,  it  is  shown  that  these  conditions  are  necessary. 

In  section  3 the  conditions  are  also  shown  to  be  sufficient. 

2.  Necessary  Conditions: 

By  the  asymptotic  stability  in  the  large  of  system  (3)  we  mean  that  the  solu- 
tions are  bounded  for  bounded  x(0) , and  that  all  solutions  x(t)  tend  to  the  origin 
as  t 4 “.  In  this  section  we  derive  restrictions  on  equation  (3)  necessary  for 
asymptotic  stability  in  the  large, 
a.  Local  Attractiveness:  The  Linear  Part 

A well-known  result  of  Lyapunov  theory  states  that  if  the  linearized  equations 
around  an  equilibrium  are  either  asymptotically  stable  or  unstable,  the  corresponding 
nonlinear  equations  are  also  asymptotically  stable  or  unstable.  Obviously  the 
linearized  equations  of  (3)  around  the  origin  are 


x ■ Ax 


(4) 


Hence,  a necessary  condition  for  asymptotic  stability  is  that  the  eigenvalues  of  the 
A matrix  have  non-positive  real  parts  and  that  its  eigenvalues  on  the  imaginary  axis 
be  simple.  Since  A is  a 2 x 2 matrix  this  only  excludes  having  two  zero  eigenvalues, 
b.  Global  Boundedness:  The  Quadratic  Part  B(x) 

If  the  stability  of  the  linear  equation  (4)  determines  the  local  properties, 
intuition  suggests  that  global  properties  of  (3)  are  determined  by  the  behavior 
of  the  quadratic  part  l.e.  B(x).  The  stability  characteristics  of  system  (2)  were 
discussed  in  detail  in  [1].  The  following  three  lemmas  derived  from  results  in  [1] 


i 


I 

l 


show  that  this  is  indeed  the  case 
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Lenrna  2.1: If  the  linear  part  of  system  (3)  Is  stable  with  A singular  and  the  quadratic 

system  Is  type  (1)  (l.e.  B(x)  i 0 at  any  point  other  than  the  origin)  then  the  system 

has  either  unbounded  solutions  or  off  origin  singularities. 

T 

Proof:  Since  A Is  singular,  let  A ■ ab  . We  consider  three  possible  cases. 

T A 

(1)  Let  B(x^)  - Xa  for  some  x^  which  does  not  lie  on  the  b^  line  (b  b^  ■ 0) . Then 
the  point 


x 


Mb^x^) 


x 


1 


Is  an  off-origin  singularity  of  equation  (3). 

2 T T 

(II)  If  B(x)  i Xa  for  any  x e p and  X e |R  let  V ■ yx  a±  where  y e |R  and  a a±  ■ 0 

• T 

V ■ yB(x)  a^  and  by  the  assumption  made  on  B(x) ,V(x)  is  > 0 for  y ■ 1 or  -1. 

Let  V(x)  > 0 for  y “ +1.  Since  V assumes  positive  values  in  the  neighborhood  of 
the  origin  the  system  is  unstable  by  Lyapunov's  first  instability  theorem. 

(III)  We  now  consider  the  case  where  B(x)  ■ Xa  only  on  bA.  Again  we  have  two 
possible  cases:  When  a “ b^,  along  the  vector  a,  Aa  - 0 and  B(a)  ■ Xa,  Hence  (3) 
behaves  as  a pure  quadratic  scalar  differential  equation  along  this  line  and  must 
be  unstable. 

When  a i b±  the  above  argument  Is  no  longer  valid  but  we  can  show  that  the 
trajectories  in  a sector  In  the  neighborhood  of  b^  must  be  directed  outwards  making 
the  origin  unstable.  A proof  of  this  is  given  In  Appendix  I. 

Lemma  2.2:  If  the  linear  part  of  system  (1)  is  stable  with  nonsingular  matrix  A 

and  the  quadratic  part  B(x)  is  of  type  (1)  then  the  overall  system  must  have  off-axis 

singularities. 

Proof:  As  was  shown  In  [1]  for  the  pure  quadratic  case  equation  (3)  can  be  written 


• 2 

xi " *1^*1 + g^xi 
• 2 

x2  - a2(v)x1  + h(v)x2 


as 
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where  v ■ — g(v)  and  h(v)  are  second  degree  polynomials  in  V 
X1 

xTGx  “ g(v)x^;  x^Hx  ■ h(v)x^  and  Ax  ■ 

v - a constant  denotes  a line  with  constant  slope  in  R*.  Assume  (without  loss  of 
generality  as  shown  in  Appendix  II)  that  *22&22  ~ al2h22  **  1116,1  the  Polynomial 

t(v)  ■ ajCvJgCv)  - a^(v)h(v)  is  cubic  and  must  have  a real  root  v^.  In  [1]  it  was 
shown  that  system  (2)  cannot  be  type  (i)  unless  h(Vg)  - 0 implies  g(v^)  i 0 and 
vice  versa. 


^(vjx 

a~7v5x~ 


t(vn) 


v v0' 

k(vn)  * 0 then  a (v_)  ■ — c-  ■ 0;  hence  a and  h have  the  common  factor 

0 20  g(vn)  2 


(v-Vq)  ®nd  Xj  = 0 along  the  line 


x0(a) 


lVj 


. ‘ai(vo) 

o g(v0)  1 " “2  “ “'“O' 

Suppose  h(vjj)  + 0 and  g(vQ)  + 0 when  t(Vjj)  » 0 then  let 


lv  0 

Then  choosing  an  - ■ \ V we  have  x,  • x,  • 0 at  x(on) 

g'v(r 


0 g(vj 


a2(v0) 

h<0  *“*  X1  • *2  ' ° “ W • 


Qualitatively,  if  A is  non-singular  this  proof  shows  that  for  some  x,  Ax  and 
B(x)  have  the  same  direction  so  that  x ■ 0 for  a suitable  choice  of  x.  The  poly- 
nomial t(v)  is  constructed  such  that  its  zero  is  the  desired  direction  v. 

On  the  basis  of  lemmas  2.1  and  2.2  it  is  clear  that  a necessary  condition  for 
global  asymptotic  stability  is  that  (3)  may  be  expressed  in  the  form 


x ■ Ax  + c xDx 


(5) 


The  following  lemma  demonstrates  that  the  system  (5)  cannot  be  globally  asymptotically 
stable  when  the  matrix  D has  real  eigenvalues.  In  [1]  it  was  shown  this  corresponds 
to  an  unstable  pure  quadratic  system  (2).  In  [2]  it  is  shown  that  for  every  quad- 

* 

Note  Oq  ■ 0 iff  A is  singular. 
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ratic  system  in  which  the  field  vanishes  on  a line  there  exists  a function 
T 

V(x)  - x Px  which  satisfies  the  Chetaev  conditions  of  instability. 

Lemma  2.3:  For  a system  given  by  equation  (5)  in  which  A is  stable,  if  D has  a 

real  eigenvalue  then  the  origin  is  not  globally  attractive. 

T 

Proof:  Let  V - x Px  where  P is  an  Indefinite  matrix.  P is  chosen  in  such  a manner 


that  if 


I>x  - (x:V(x)  > 0} 


T T 

then  c L ^ where  L 2 « (x:x  (PW-D  P)x  > 0).  Further  c^  l L Thus  a closed  region 
V > 0 is  embedded  in  an  open  region  V > 0 and  the  system  is  unstable.  [For  a more 
detailed  proof  the  reader  is  referred  to  [2].] 

c.  Other  Necessary  Conditions:  The  three  previous  sets  of  conditions  in  this  sec- 
tion have  had  rather  obvious  intuitive  motivation.  Asymptotic  stability  of  the 
linear  part  assures  (local)  asymptotic  stability  of  the  nonlinear  system  as  well. 
Stability  of  the  quadratic  part  is  required  to  assure  the  boundedness  of  all  solu- 
tions with  bounded  initial  conditions.  The  final  conditions  are  related  to  both 
linear  and  quadratic  parts.  They  are  needed  to  ensure  the  absence  of  off-origin 
singularities  and  a well  behaved  interaction  between  the  linear  and  quadratic  parts 
far  away  from  the  origin. 

Lemma  2.4:  The  system  described  by  the  equation  (5)  with  A asymptotically  stable  has 
off-origin  singularities  if  and  only  if  the  matrix  A has  real  eigenvalues. 

Proof : Equation  (3)  can  also  be  expressed  as 

x - A[I  + «TxA  *D]x 

since  A is  nonsingular.  If  A has  an  eigenvalue  A and  corresponding  eigenvector  x. 


Unless  A D has  an  eigenvalue  of  multiplicity  2,  whose  associated  eigenvector  is 
in  the  direction. 
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- is  an  off-origin  singularity. 

-X(c  *x) 

T 

If  c Xj  • 0 we  choose  the  eigenvector  corresponding  to  the  other  eigenvalue. 

If  A has  only  complex  eigenvalues  then  no  vector  x^  exists  such  that 
T -1 

c xA  Dx  ■ -x. 


f 


Since  A is  nonsingular  this  implies  that  there  are  no  off-axis  singularities. 

In  addition  to  precluding  off-origin  singularities  the  condition  that  A ^D 
have  complex  eigenvalues  affects  the  Interaction  of  the  linear  and  quadratic  parts 
of  (5)  far  away  from  the  origin.  Since  the  quadratic  part  B(x)  ■ 0 along  the  line 
c^ , the  influence  of  the  linear  part  is  felt  in  the  neighborhood  of  this  line  even 
as  x -*■«».  The  following  lemma  addresses  the  critical  interaction  of  linear  and 
quadratic  parts  near  c^. 

Lemma  2.5:  If  the  matrices  A and  D have  complex  conjugate  eigenvalues,  then 

2 

for  every  vector  x £ R there  exists  a X such  that  x is  an  eigenvector  of  the  matrix 
(A+XD) . 

Proof : Since  yx  + A ^Dx  for  all  y e R we  know  that  yAx  i Dx  for  all  y e R and  all 

2 2 
x £ 1R  . Hence  Ax  and  Dx  are  linearly  Independent  and  span  (R  and  x can  be  expressed 

as  a linear  combination  of  the  two. 

a Ax  + BDx  ■ x 


or 


[A  + 8/a  D]x 


1 

a 


x 


Hence  choosing  X - e/a  the  result  follows. 

Since  any  direction  is  an  eigenvector  for  a suitable  linear  combination  of  A 
and  D we  have,  in  particular 


[A  + x0dJca  ■ Yc± 


W 


Unless  A 1D  has  an  eigenvalue  of  multiplicity  2,  whose  associated 
in  the  c^  direction. 


eigenvector  is 


i 


i 


I 


r 


* 

I 


where  is  the  unit  vector  which  is  the  orthogonal  complement  of  c.  Then  writing 


(5)  as 


where 


x - [A  + 1qD]x  + o(x)Dx 


o(x)  « c x - X, 


the  system  (6)  behaves  essentially  like  a linear  system 


x - [A  + XqD]x  (7) 

at  points  far  from  the  origin  and  In  the  neighborhood  of  the  c direction.  In 
1 T 

particular  if  x = — c + e c then  o(x)  *■  e c c - X is  a constant  and  the  field 

E1  1 2 2 Ci 

of  system  (6)  tends  towards  that  of  (7)  as  0. 

Equations  (6)  and  (7)  provide  the  analytic  expression  of  our  understanding  that 
system  (5)  Is  approximately  linear  not  merely  in  the  neighborhood  of  the  origin  but 
also  for  large  values  of  x on  the  direction  where  the  quadratic  part  vanishes. 

The  obvious  question  then  is  whether  the  linear  system  (7)  is  stable.  This  is 
answered  by  considering  the  eigenvalue  y corresponding  to  the  eigenvector  c of  the 
matrix  [A  + X^  D] . The  condition  y < 0 is  necessary  for  global  asymptotic  stability 
is  demonstrated  by  the  following  lemma. 

Lemma  2.6:  If  system  (7)  has  an  eigenvalue  Y > 0 corresponding  to  the  eigenvector 
c^  then  system  (6)  has  unbounded  solutions. 

Proof:  We  will  construct  a region  which  is  disconnected  from  the  origin  and  which  is 
invariant:  i.e.  any  trajectories  starting  on  its  interior  must  remain  so  for  all 
time.  By  this  construction  we  will  have  violated  the  necessary  condition  that  a 
globally  asymptotically  stable  equilibrium  state  be  the  only  proper  Invariant  set 
of  Bt^. 


■ 


Consider  system  (6)  when  y > 0 and  assume,  with  no  loss  of  generality,  that 
XQ  < 0.  Choose  a point,  x^,  such  that  ||Xjll  » 1 and  o(x^)  - 0 where 


-8- 


x - — c.  + e0c,  and  then  define  x_  - c - e c.  The  region  in  question  will 

1 x z i.  x z 

be  the  triangular  wedge  bounded  by  the  lines  {Xx^|x  > 0}  and  {x^  + XX2IX  > 0}  as 
shown  in  Figure  1. 

To  show  invariance, 
it  suffices  to  show 
that  the  field  on 
the  entire  boundary 
is  directed  towards 
the  interior  of  the 
region. 

At  the  point 
x^  the  system  looks 
like  equation  (7) 

\ 1 

o(x)-0 

[A  - X^Jxj.  - Y l*1  cx  + e3c  . 

For  and  sufficiently  small  the  field  is  directed  toward  the  interior  of  the 
region.  For  large  values  of  X along  either  bourding  line,  o(x)  becomes  very  large 

and  the  field  is  governed  by  o(x)Dx.  As  established  in  [1]  one  half  of  the  line 
through  the  origin  in  the  direction  is  an  attractive  set  and  without  loss  of 
generality  we  can  assume  that  the  constructed  region  contains  part  of  this  attractive 
set.  Since  the  pure  quadratic  field  is  directed  towards  the  c^  line  the  field  in 
(6)  is  the  resultant  of  the  radially  outward  linear  part  and  the  tangentially 
Inward  quadratic  part  and  is  obviously  directed  into  the  interior  of  the  region. 

3.  Sufficient  Conditions; 

In  the  previous  section,  it  was  shown  that  conditions  (l)-(lv)  are  necessary 

1 

for  the  asymptotic  stability  in  the  large  of  equation  (3) . In  this  section,  it 


is  shown  that  these  conditions  are  also  sufficient. 


Lemma  3.1:  If  in  equation  (6),  c^  is  an  eigenvector  corresponding  to  a negative 
eigenvalue  of  the  matrix  A + AQD,  all  solutions  of  the  differential  equation  (3) 
are  bounded. 

Proof : This  proof  is  based  upon  the  behavior  of  the  solution  in  three  distinct 
2 

regions  in  (R  : a region  where  the  quadratic  part  of  (6)  dominates  the  field, 

which  we  denote  N;  a region  where  the  linear  part  of  (6)  dominates,  denoted  L; 

and  the  interface  between  them,  denoted  T.  These  regions  will  be  defined  with 

reference  to  three  lines: 

A 

b * c + e c^  (where  e > 0 is  any  arbitrarily  small 


d **  c - e c 


o(x)  » c x - = 0 


constant  of  order  — ) 
*0 


The  vectors  b and  d define  a conical  region  containing  the  l*ne  c^  and  a(x)  » 0 
is  an  affine  line  parallel  to  c^.  Specifically,  we  will  define 


N - {x  £ IR2 


xTbdTx  2 0 and  \\  x l\  5 Rq} 


as  the  portion  of  excluding  a disc  of  radius  and  the  cone  containing  the 
line  c. . 


L - (x  e |R2 


|o(x)|  ||  D l|  < Y1  < Y.  l\xtl  > Rq> 


is  an  open  strip  around  the  line  o(x)  ■ 0 upon  which  the  linear  field  (7)  dominates 
in  (6) . 


A / _2 

T - (x  e ft 


xTbdTx  < 0,  l|xl\  > Rq) 


defines  a conical  region  around  the  cx  line  and  contains  L for  sufficiently  large 
Rq.  The  sets  are  depicted  in  Figure  2. 

Any  point  x e L lies  close  to  the  direction  so  that  the  field  may  be 
arbitrarily  approximated  (as  Rq  Increases)  by 
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f 

. 

I 

. 


(A+XqD)x  + o(x)Dxa#  -yx  + jj^jj  Dx 

which  Is  directed  radially  Inwards.  Hence,  no  solution  which  remains  In  L can  be 
unbounded. 

If  bL+,  d±+  and  b±  and  d±  denote  the  half  line  boundaries  of  the  portion 

of  T containing  the  attractive  and  repulsive  halves  of  the  c^  line  respectively 

(see  Figure  2)  then  the  field  Is  toward  the  Interior  of  N along  b^  and  d^  and 

toward  the  interior  of  T along  b^  and  d^+.  Along  all  trajectories  in  N the  field 

T 

is  also  governed  by  the  quadratic  part  c xDx  with  a small  perturbation  due  to  the 

linear  part  Ax.  Defining  tan  6(x)  - J** — * — — < .a  (where  a - Max||Ax|j, 

T IcTxDx«  6|cTx|  6c1R0  |jx|/-l 

6 » Min  Dx  and  |c  x|  > c Rq  In  N)  we  can  bound  the  maximum  angular  deviation  of 
Rxfl 

the  perturbed  trajectory  from  that  of  the  pure  quadratic  solution.  Since  the  arc 
length  of  the  unperturbed  trajectory  Is  proportional  to  the  maximum  deviation 
of  the  perturbed  trajectory  is  bound" j by 


I - — s—  ds  - ( <*here  T is  the  path  of  the  trajectory  in  N) 

J 6ClRQ  0 

where  a ^ is  a constant.  The  unperturbed  trajectory  leaves  N at  a finite  point 
along  b + and  hence  the  perturbed  solution  must  do  so  as  well.  Hence  no  solution 
in  N can  be  unbounded. 

Finally  in  T,  the  quadratic  part  predominates  near  the  boundary  and  weakens 
as  a(k)  -*-0  - i.e.  as  L is  approached.  However,  near  L,  where  the  effect  of  the 
quadratic  term  is  weak,  the  linear  field  is  directed  radially  inward  as  demonstrated 
before.  Thus,  on  the  interior  of  T the  field  is  the  resultant  of  a tangential 
component  and  a radially  directed  component.  Hence,  no  solution  can  become  unbounded 
within  T. 

We  have  shown  that  unbounded  solutions  cannot  occur  within  the  set  Tt;LuN.  Since 


J 


the  set  {x|  If  x ||  < RQ}  which  is  its  complement  is  itself  bounded  this  suffices  to 
preclude  the  possibility  of  any  unbounded  solutions. 


I 


f 


Since  all  trajectories  are  bounded,  the  solutions  oust  approach  an  invariant 

set.  Under  conditions  (i)-(iv)  the  system  has  an  attractive  origin  but  no  other 

equilibrium  states  and  cannot  sustain  any  limit  cycles.  Hence,  the  origin  is  the 

2 

unique  proper  Invariant  subset  in  ft  and  all  solutions  tend  to  zero.  Hence,  the 
system  is  globally  asymptotically  stable. 

Conclusion: 

This  is  a sequel  to  the  earlier  report  on  pure  quadratic  systems 
considers  the  complexities  Introduced  by  the  addition  of  a linear  term.  The 
diversity  of  qualitative  uchavior  evinced  by  such  systems  has  led  us  to  limit 
this  investigation  to  a specific  question  - global  asymptotic  stability.  Necessary 
and  sufficient  conditions  for  such  stability  are  derived  here  taking  advantage  of 
the  relatively  simple  structure  of  the  pure  quadratic  systems  discussed  in 
These  conditions  may  find  application  in  such  areas  of  systems  theory  as  adaptive 


control  and  bilinear  systems. 
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• T 

Fact:  If  x ■ ab  x + B(x)  where  B(x)  ■ xa  if  and  only  if  x ■ then  the  origin  la 
unstable. 

Proof : In  the  body  of  the  paper  it  was  shown  that  the  origin  Is  unstable  when 

a ■ b^ . Hence,  we  assume  that  a and 

are  linearly  independent.  It  can  be 

shown  that  on  one  side  of  the  line 

b±  there  is  a curve  C (asymptotically 

approaching  bx  as  tyx\|  decreases)  on 

which  the  component  of  B(x)  in  the 

direction  of  the  vector  'a'  is 

T 

exactly  equal  to  (-b  x)a.  Thus 
the  resultant  vector  must  be  in 
the  direction  b^  or  -bx  at  every 
point  on  the  entire  curve  C. 

If  vectors  b^  and  b^  are  defined 
as  ~ 

A 

\>1  - bx  + e b 
. A 

b2  " “bi  + E b 

the  curve  C passes  through  every 
triangular  wedge  bounded  by  b^ 
and  bx  and  b 2 and  -bx  as  e decreases. 

Further,  in  the  neighborhood  of  the  Figure  3 

origin  all  trajectories  on  b1»b2,bi  and  ~bi  are  directed  towards  the  interior  of 
the  corresponding  triangular  wedge.  Hence,  if  the  component  of  B(x)  in  the  bA 
direction  is  radially  outward  on  one  of  the  triangular  wedges,  the  origin  must 
be  unstable. 


I 
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It  remains  to  show  that  If  the  resultant  is  radially  Inward  on  one  side  It 
must  be  radially  outward  on  the  other.  To  do  this,  parametrize  the  curve  C by 
c(w)  where  c(0)  ■ 0.  Since 


11m 


■?.W 


and 


lim  -c(-oi)  m i 
\\  c(-w)»  \\  b^ 


and  B(Xc(u>))  * X B(c(a>))  we  may  study  the  properties  of  B(x)  over  an  arbitrarily 
small  neighborhood  of  b^.  t 


Consider  the  mapping  v»-> 


h(v) 

g(v) 


A X2 


used  in  reference  (1).  v ■ — ; 


1 T 1 T 

g(v)  - — 2 x Gx5  b(v)  - — j x Hx.  In  the  vicinity  of  Vq  (slope  of  b^)  there  is 

change  in  the  sign  of  ihe  tangent  to  that  curve  at  h(Vg)/g(Vg).  This  implies 


no 


(ref.  [2])  that  there  exists  another  direction  v^  for  which 

h(vx)  h(vQ) 

g(Vj)  g(vQ) 


But  this  contradicts  the  condition  that  B(x)  - Xa  only  on  b 


Hence  must 

g(v) 


take  on  a local  extremum  in  the  neighborhood  of  v^  in  which  case 


b±  B(c(to))  < 0 


e > to  > 0 


and 


-b  B(c(-u»))  > 0 


so  that  if  the  resultant  is  radially  inward  on  one  side  it  is  radially  outward  on 
the  other  and  the  system  is  unstable. 
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Appendlx  II 

We  will  show  that  if  x - Ax  + B(x)  and 


*22822  “ ®12h22  " 0 


then  the  system  has  an  off-origin  equilibrium. 


Since 


12  _ ^22 

*22  h22 


the  vector 


12 

*22 


and 


g22 

h22 


are  in  the  same  direction  in  R . But 


"o 

” al2~ 

/ 

— — 

0 

\ 

g22 

A 

1 

and  B 1 
\ 

1 

L j 

CM 

CM 

<0 

mm 

/ 

h22 

and  hence 


- XB 


we  Jiave 


x - 0 at  the  point 


and  choosing  y ■ -1/X 


’ 


mm 


